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Method and apparatus for image processing. 
The invention relates to a method and apparatus for processing images. More 
specifically, the invention is concerned with detecting and identifying specific 
features or types of features in images, and classifying objects and regions in 
5 images. 

The invention is especially useful for the detection of man-made structures in 
images, which has applications, for example, in photo-interpretation, 
cartography and surveillance. 

Methods of identifying features such as lines and edges in images are known. 

10 Simple methods of linking edges into object boundaries are also known, but 
these are generally successful only for relatively simple scenes. It is also 
known to use model-based techniques for identifying specific features in 
images, but these techniques require prepared reference models, which may 
not necessarily be available for all features. Other approaches attempt to infer 

15 3-D structure, but are often not reliable and may require multiple views of a 
scene. 

Furthermore, it is especially difficult to detect man-made structures and urban 
areas in images, for a number of reasons. For example, contrast between 
buildings and surrounding structures can be very low, causing low-level 

20 segmentation and edge extraction to be fragmented. In addition, small 
structures, such as cars, trees or rocks adjacent to large structures, such as 
buildings, can cause further problems at a segmentation stage and give rise to 
noisy and perturbated, or phantom boundaries. The problem is often further 
complicated by the existence of shadows and markings. 

25 In the specification, references to an image include a whole image and a 
region of an image. 

Aspects of the invention are set out in the accompanying claims. 
According to a first aspect of the invention, there is provided a method of 
analysing an image comprising performing a Hough transform on points in an 
30 image space to an n-dimensional Hough space, selecting points in the Hough 
space representing features in the image space, and analysing m of the n 
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variables for the selected points, where m is less than n, to derive information 
about the features in the image space. 

More specifically, feature points in an image are extracted, for example, using 
feature detecting means such as an edge or corner detector, or other similar 
5 feature detecting means. The feature points undergo the Hough transform, 
and the image of the feature points in the Hough space are analysed. 
As a result of the invention it is possible to draw more reliable conclusions 
about features in images. The technique does not require a model of a 
structure and handles boundary and segmentation problems well. 
10 Embodiments of the invention will be described with reference to the 
accompanying drawings of which: 

Fig. 1 is a flow diagram illustrating a method according to an embodiment of 
the invention; 

Fig. 2 is a diagram illustrating the Hough transform; 
15 Fig. 3 is another diagram illustrating the Hough transform; 
Fig. 4 is an example image showing man-made structures; 
Fig. 5 is an image showing the results of an edge detection process on the 
image of Fig. 4; 

Fig. 6 is an image showing the result of a Hough transform performed on the 
20 points in the image of Fig. 5; 

Fig. 7 is an image showing the results of an inverse Hough transform based on 
Fig. 6; 

Fig. 8 is a graph showing the results of a mapping from Hough space based on 
Fig. 6; 

25 Fig. 9 is a schematic diagram of an apparatus according to an embodiment of 
the invention; 

Fig. 10 is a flow chart illustrating a second embodiment of the invention. 
The first embodiment will be described with reference to the example image 
region shown in Fig. 4, which includes a number of buildings, ie man-made 
30 structures, and the flow diagram in Fig. 1. 
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In this embodiment, the method is performed using a computer, illustrated 
generally by the reference number 200 in Fig. 9, programmed to execute the 
method using suitable software. The computer 200 includes a processor 210, 
image display means in the form of a display 220 and image input means in 
5 the form of a scanner 230. 

As a first step, step 50, a digital image is input to the computer. In the present 
embodiment, the image is input using a scanner 230 scanning a photograph, 
although other input means may be used, or an image may be retrieved from 
data storage means. 

10 At step 100, the image is divided into sub-regions. Here, the regions are 
blocks of 256x256 pixels. The blocks may be overlapping or non- 
overlapping. In this embodiment, the blocks are non-overlapping and 
adjoining. Fig. 4 shows one block. 

The next step, step 1 10, is the detection of feature points in each block. In the 
15 present embodiment, the positions of feature points of potential interest are 
identified in the image using an edge detector, which identifies discontinuities 
in pixel values. Various edge detectors are known and can be used; in this 
embodiment a Canny edge detector is used. 

The edge detection step classifies each pixel in the image regions as either a 
20 one of feature points or not a feature point and produces a binary image 
accordingly. Fig. 5. is a binary image showing feature points in the image 
region of Fig. 4. 

In step 130, the feature points in Fig. 5 are represented in Hough space using a 
Hough transform for line detection. 

25 A general discussion of the Hough transform follows. 

Figs. 2 and 3 illustrate the Hough transform. Hough transform is used to 
detect curves that can be expressed in a parametric form. The Hough 
transform projects feature points in image space (Fig. 2) to points in 
parameter space also called Hough space (Fig. 3). Each point in the image 

30 space "votes" for all sets of parameters (e.g. points in the hough space) that 
generate parametric curves passing through that point. Referring to Fig. 2, for 
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a point A(x,y), consider one line passing through that point. The 
perpendicular from the line to the origin has a length p and is at an angle 9 
from the x axis, and the variables p and 8 define that line in Hough space, as 
shown by the dot in Fig. 3. The Hough transform of a point in image space 
5 takes all the possible lines passing through a point A(x,y) and plots the 
resulting (p,0) parameter values in Hough space. 

Fig. 3 shows the boundaries in Hough space. The horizontal axis corresponds 
to the 0 values which vary from 0 to 2tc, as can be understood from Fig. 2. 
The vertical axis corresponds to the p values and varies between a minimum 
10 of a and a^Jl in cycles for each nil radians, where 2a is the width of the 
block shown in Fig. 2. 

In the Hough transform, each feature point in image space is transformed to 
points in Hough space. Processing of plurality of point in image space results 
in accumulation of points at some p, 0 locations in the Hough space, and co- 

15 linear points in the image space result in peaks, or local maxima, in the Hough 
space. The local maxima are extracted and then an inverse Hough transform 
is applied to map back into the image space to derive each line corresponding 
to each respective local maximum. The resulting lines are of infinite length , 
and so the lines are combined to get an idea of the shapes corresponding to the 

20 feature points in the original image space. 

In this specification, the dimension of the Hough space relates to the number 
of parameters defining the parametric curve. Thus, in the case of a Hough 
transform for detecting lines, a point (x, y) in image space is mapped onto a 
point (p, 0) in two-dimensional Hough space. 

25 Fig. 6 illustrates the result of the Hough transform performed on the feature 
points in Fig. 5. The intensity values in Fig. 6 correspond to the accumulated 
values so that peaks in Hough space are shown as the brightest points (local 
peaks being the locally brightest points). 

In step 140, landmark points are selected. First, the local maxima of the 
30 accumulator values in the Hough space are detected. Fig. 7 illustrates the 
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lines in image space corresponding to the local peaks in Hough space. The 
accumulator values for the local maxima are compared with a threshold. In 
this embodiment, the threshold is a predetermined percentage (80%) of the 
global maximum accumulator value, for the block. The (p, 9) values for the 
5 local maxima which exceed the threshold are selected as a subset, which is 
called the landmark subset. The threshold may vary. For example, other 
percentages can be used, and the global maximum may instead be for the 
whole image, or the threshold may be determined in other ways. 
The local maxima in the landmark subset are then processed further, in step 

10 150. In this embodiment, the (p, 0) values in the landmark subset are 
projected onto the 0 axis and accumulated. This results in a graph as shown 
in Fig. 8 for the present example. The result of the projection is analysed to 
determine characteristics about the lines corresponding to the landmark 
subset. More specifically, in this embodiment, the relative angles of the lines 

15 for the landmark subset are analysed. As can be seen from Fig. 8, the highest 
peaks correspond to prominent directions of 46 and 136°. This indicates a 
large number of lines intersecting at 90°, which suggests a man-made 
environment such as urban area. Similarly, a large single peak suggests 
parallel lines which generally occur in a man-made environment. 

20 Other projections onto other variables could be used. For example, the 
projection could be onto the p axis. 

A second embodiment will now be described with reference to the flow chart 
in Fig. 10. The flow chart in Fig. 10 is the same as the flow chart in Fig. 1, 
except that there is an additional step, step 120, between steps 110 and 140. 
25 The second embodiment is similar to the first embodiment, apart from 
differences discussed below. 

In the step of detection of feature points, step 110, the edge detection in this 
embodiment assigns values indicating whether a feature is present at given 
location, to points in the image. For example, an edge confidence measure 
30 taking values between 0 and 1 can be used, higher values indicating higher 
confidence in existence of an edge, resulting in a probability map for edge 
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feature points. Other features can be also used, for example corners or colour 
features. 

In step 120, after step 110, a statistical analysis of the feature points detected 
in step 110 is performed for each block. More specifically, the confidence 
5 values are made into a histogram for each block. Also, the density of feature 
points for each block is calculated. Other statistical information relating to 
the distribution of feature points in the blocks may also be derived. This 
statistical information is used in subsequent processing, in selecting landmark 
points. 

10 A threshold for detecting landmark points can be derived as follows. A large 
number of randomly selected representative configurations of feature points 
are generated for a given block size. For each configuration, the Hough 
transform is performed, and a histogram of accumulated values in Hough 
space is determined. The histograms for all the configurations are combined 

15 (eg averaged) and used to determine a suitable threshold for detecting 
landmark points. 

For example, a value equivalent to the highest value observed in the combined 
histogram may be used as this threshold. 

The combined histogram is a reflection of the image in Hough space arising 
20 from general random image behaviour, in the particular application, eg for the 
specific block size. In other words, this approach gives an indication of peaks 
in Hough space that occur randomly and do not necessarily correspond to 
landmark points in an image, ie noise. The threshold can be varied according 
to the level of confidence required in the landmark point extraction. For 
25 example, a higher threshold, such as three times the highest value in the 
combined histogram, increases the exclusion of landmark points that occur by 
chance, but on the other hand it increases the risk of excluding possible points 
of interest. Similarly, a reduction of the threshold, such as to 80% of the 
highest value in the combined histogram, increases the likelihood of returning 
30 points of interest correctly, but also increases the risk of returning false points. 



7 



Other statistical information as described above can be used in deriving the 
threshold information. For example, the randomly selected representative 
configurations of feature points may take into account the density of feature 
points for the image region under consideration, that is, they have a similar 
5 density, and similarly for the other statistical information as appropriate. 

In order to simplify processing, the modelling for generating reference 
histograms from the randomly generated images and corresponding thresholds 
can be performed offline and the histograms or threshold values stored, for 
example, in a look-up table. 

10 The embodiments described above use a Hough transform for detecting lines. 
The Hough transform can also be used to detect other regular curves such as 
circles, ellipses, and even in applications where a simple analytic description 
of a feature is not possible, although in such cases the computational 
complexity increases. The invention can also be used in these other 

15 applications of the Hough transform principle. For example, the invention 
could be used in analysing the relationship between the centres of circles, 
where a regular spacing of circle centres suggests human influence. 
The invention can be implemented by any suitable apparatus, such as a 
suitably programmed computer which is either a dedicated application 

20 specific computer or a computer supplied with suitable software to execute 
the method or implemented in hardware such as a specially designed 
integrated circuit. The apparatus may include image input means such as a 
scanner, digital camera, internet link, or a connection to image storage means 
such as a data carrier. 



